Randomly dilute Ising model: A nonperturbative approach 
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The N- vector cubic model relevant, among others, to the physics of the randomly dilute Ising 
model is analyzed in arbitrary dimension by means of an exact renormalization-group equation. 
This study provides a unified picture of its critical physics between two and four dimensions. We 
give the critical exponents for the three-dimensional randomly dilute Ising model which are in good 
agreement with experimental and numerical data. The relevance of the cubic anisotropy in the 
O(N) model is also treated. 
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Understanding the effects of a weak quenched disorder 
on magnetic systems is one of the great challenges of sta- 
tistical physics. In this context, much attention has been 
paid to the Ising model with randomly distributed non- 
magnetic-impurities, i.e. the random Ising (RI) model 
(seeEl anda for reviews) described by the following Hamil- 
tonian : 



H — — J EiEjSiSj 
(id) 



(1) 



with J > 0. In Eq. ([!]), the Si are Ising variables and the 
Si are quenched, uncorrelated, random variables taking 
the value 1 with probability p and with probability 
1 — p. According to the Harris criterionH, the critical 
behavior of a pure system is altered by the disorder if its 
specific-heat exponent a is positive, as it is the case for 
the three-dimensional Ising model for which a ~ G1109, 
(3 ~ 0.326, 7 ~ 1.24, v ~ 0.630 and r) ~ 0.0335 (seJ). 

Experimentally, there now exist several compounds 
corresponding to a RI model. They essentially con- 
sist in crystalline difluoride of iron or manganese where 
FeF2 or MnF2 is replaced by a nonmagnetic material 
such as ZnF2. Experiments performed on Fe p Zni_ p F2 or 
Mn p Zni_ p F2 with different values of the concentration 
p have revealed scaling laws with exponents clearly dis- 
tinct from those of the pure Ising model. One has (seea), 
for instance, for Fe p Zni_ p F 2 : —0.12 < a < —0.06, 
0.34 < /3 < 0.38, 1.28 < 7 < 1.4 and 0.68 < v < 0.72. 

Numerically, early Monte Carlo simulations had con- 
firmed the relevance of disorder for the Ising model (seeo) . 
The situation was however unclear since the critical ex- 
ponents seemed to vary with p. Recently, Ballesteros et 
al.tj, by a careful infinite volume extrapolation of the data 
based upon an analysis of the correction-to-scaling terms, 
have reached the conclusion of dilution-independent crit- 
ical exponents : a = -0.051(16), = 0.3546(28), 
7 = 1.342(10), v = 0.6837(53) and rj = 0.0374(45) for 
0.4 < p < 0.9. 

From a theoretical point of view, the continuous field 
theory relevant to the study of the RI model is obtained 
by considering N replica of the model (Q) and by taking 



the average over all possible realizations of the Gaus- 
sian distributed disorder. This leads to the Ginzburg- 
Landau- Wilson action for the TV-vector model with cubic 
anisotropytffi : 



S = J d d xKJ2[(d^ a ) 



2 ' m^ a ] 



(2) 
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with u\ oc — (1 — p) < and 1*2 > 0. The physics of the 
RI model is then obtained through the N — > limit of 
the renormalization group (RG) equations obtained from 
action (Q). The determination of these RG equations 
has been the subject of an intense activity, mainly per- 
formed within perturbative frameworks. Different renor- 
malization schwajiflsJiaive been used : minim al . aubtoaflti ea 
(MS) schemeQBElll3'til or massive scheme^HlljIlalla'tZl 
directly in three dimensionSr-where the (3 functipns have 
been computed up to fivetil and six loopsla respec- 
tively. In each scheme, two different kinds of expan- 
sions for the corresponding (3 functions have 
sidered : i) e-expansion (actually -^e-expansior 
in the MS scheme and pseudo--y/e-expansiont3 in the 
massive scheme it) expapsion iip-the coupling constants 
in three dimensionsO^Ea£3ll3'0. It has, however, ap- 
peared that the series obtained in case i) are non-Borel 
summable and thus inappropriate for the— mmntitative 
evaluation of critical exponents, b^haviomE3Ejl, a result 
anticipated by several authorsE3Ell within an .analysis 
of the zero-dimensional problem (see howeverE3). In 
case it), generalizations of Pade-Borel resummation tech- 
niques for several variables have been considered. In 
the MS scheme. Chisholm-Borel resummation scheme has 
been usedBEjCJ. This method has led to satisfying results 
at low order but has also appeared to be unstable with 
respect to the order of approximation. Indeed, the five- 
loop calculation does not allow to obtain real coordinates 
for the RI fixed point. This point has been wj^^ed to 
the possible Borel-non-summability of the series3tL3'El In 
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the massive scheme, Pelissetto and VicariEl, using a pro- 
cedure of resummation formulated within .-a, subsequent 
analysis of the zero-dimensional problerr£3, have em- 
ployed Pade-Borel and conformal Pade-Borel resumma- 
tion techniques within a six-loop RG calculation. They 
found a fixed point with exponents 7 = 1.330(17), v = 
0.678(10), r] — 0.030(3) that are in rathej^imod agree- 
ment with previous perturbative result s! 1 ^! 14 !] 16 ! and with 
the present-day numerical data. Let us however mention 
that, even in this case, the location of the fixed point 
is found to be calculation-dependent. This has led Pelis- 
setto et al. to consider the possibility of an optimal trun- 
cation, possibly shorter than that considered ir£3, a fact 
that could constitute a limitation of the method. 

Despite the good agreement between experimental and 
numerical data and the perturbative results, we would 
like to stress upon the need for a complementary ap- 
proach, not relying on perturbative expansions and this, 
for several reasons. The first one is that the perturba- 
tive field theoretical treatment implies the use of sophis- 
ticated resummation techniques which, if they proved 
to lead to safe results in the case of pure magnets, are 
far frorftihftiijig; .under control in the case of disordered 



systemsl 



The second reason for the use of a non- 



perturbative treatment is that there exist cases where 
the perturbative approach apparently seems to provide 
incorrect results. For instance, in the context of frus- 
trated magnets, six-loop perturbativp-RG computations 
predict the existence of a fixed pointEa and, thus, a uni- 
versal behavior for their critical physics. The experimen- 
tal data, although displaying scaling laws, are rather in 
favor of a nonuniversal behavior, the critical exponents 
varying from one material to the other. Computations 
based on the exact RG formalisrruhpfe shown to be able 
to reproduce this phenomenologyE§t3. Finally, there are 
situations where only specific techniques work well. For 
example, in two dimensions, one can use integrability or 
conformal invariance which, however, only operate in this 
dimension. All these features plead toward the use of a 
method that does not rely on the treatment of divergent 
series, provides a systematic way to improve the results 
and is able to tackle with the physics in any dimension. 

A Wilson-likeEZI nonperturbative treatment of (Q) has 
been first performed by Newman and RiedelEa on the ba- 
sis on the scaling-field (SF) method. However, due to 
the complexity of its practical implementation, this ap- 
proach has been given up and perturbative approaches 
largely preferred. We used here a nonperturbative ap- 
proach, .based on the concept of effective average action 
rfc[0]E3H which appears as a very efficient tool. The 
quantity Tk [</>] has the meaning of a coarse-grained free 
energy at the scale k in the sense that it only includes 
the fluctuations with momenta q > k. It is a function of 
the effective order parameter 4> at the scale k, the ana- 
log of a magnetization at this scale. In the low-energy 
limit k — > 0, all modes are integrated out so that Tk=o 



identifies with the usual effective action T. At the mi- 
croscopical scale A, no fluctuation has been integrated 
and Tk=A identifies with the microscopic action S. The 
fc- dependence of Tk is controlled by an exact evolution 
equatiorH : 



dT 



(ri 2) 



■Rk) 



dt 



(3) 



where t — ln(fc/A). The trace has to be understood as 
a momentum integral as well as a summation over inter- 
nal indices. In Eq.Q, T^\ the second derivative of Tk, 
is the exact field- dependent inverse propagator and Rk is 
the running infrared cutoff that suppresses the propaga- 
tion of modeSj-ss/ith momenta q < k. A convenient cutoff 
is provided byo : Rk(q) — Zq 2 /[enp(q 2 /k 2 ) — 1], where Z 
is the field renormalization. The effective average action 
Tk is a functional invariant under the symmetry group of 
the system and thus includes all powers of all invariants 
- and their derivatives - built out from the average order 
parameter. Thus, Eq.^) is a nonlinear functional equa- 
tion, too difficult to be solved exactly so that, in practice, 
one has to deal with a truncated form of Tk . We choose 
here to perform an expansion of Tk around its minimum, 
keeping a finite number of monomials in the invariants 
and including the most relevant derivative termsS3. We 
first consider the following truncation : 



Tk = f d d xl^2^(d^ a ) 2 + Ul 

N 1 
o=l J 



(4) 



where {ui,U2,K, Z} are the fc-dependent coupling con- 
stants that parametrize the evolution of the model with 
the scale k. The action (|J) displays a discrete Zjv per- 
mutation symmetry. For Ami + u 2 > 0, the minimum 
of action M) is given by a replica symmetric configura- 
tion defined by 0™ m = \[k for all a. Expression (^) 
thus corresponds to a quartic expansion of Tk around 
the minimum 0™ m . There is no residual symmetry in 
the minimum so that the symmetry breaking scheme is 
given by a fully broken TL^ group. The spectrum of 

the theory in the broken phase is given by the eigen- 

(2) 

values of T k at the minimum and consists in a singlet 
of mass m 2 — 8k(Nu\ + U2) and a (N — 1) uplet of mass 
777,2 — 8K7i2- When 7i2 = 0, this (N — l) uplet corresponds 
to the Goldstone modes of the 0(A)-invariant model that 
gets here a mass proportionnal to the 0(A)-symmetry- 
breaking-coupling-constant u 2 ■ 

Using Eq.(g) and Eq.@, we have derived the flow 
equations for the dimensionless renormalized coupling 
constants n r , u\ r , u 2r , and Z : 
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d t K r 

d t Uir 
d t U 2r 



— (d— 2 + 77) /c r — 
(d - 4 + 2i])u lr - 



2v d 



N(Nui r + u 2r ) 
2v d 

N 2 K r 



3(Nu lr + u 2r )L d 1 [m( r } + (N - l)(Nu lr + 3u 2r )L d [m 2 ri 
- Nui r L([ml r ] + Nu lr Lf[ml r ] + 2>&n r (Nu lr + u 2r f L d 2 \m\ r \ 



+ 4(7V- l)N 2 K r u( r + 24N(N - l)K r u lr u 2r + 36n r u 2 r )L d [m 2 2r ] - 8K r {Nu lr + 3u 2r ] ) 2 L d A [m\ r ,m 



(5) 



(d - 4 + 2r?)u 2 r 



2vd_ 



Nu lr L d [rn{ r ] - Nu lr L d [m 2 2l ] + 36(A - 2)n r u 2 2r L' 2 l [m 2 



+8fi r (Nu lr + 3u 2 r) 2 i?.l[ m lr> m 2r 

256n r Vd 



-d t In Z 



Nd 



(TV - l)u 2 2r M d 2 [Q, m 2 2r ] + (Nu lr + U2r ) 2 M 2 d 2 [0, m? r ] 

I 



where l/vj 



\mj r , 



2 d+1 n d / 2 r(d/2). 



In 



L d Jm 2 



and M n u n 2 



[m\ r ,m 2 2r ] 



Eqs.©, 



ril,n2 L 

the "threshold functions"^ that encode the nonpertur- 
bative aspects of the RG equations @. We show how 
this set of equations provides a qualitative and, partly, 
quantitative description of the physics of model (g) 
between d = 2 and d = 4. 

Around d = 4, the fixed point coupling constants u lr 
and ii 2r are expected to be of order e = 4 — d. One 
can thus obtain simplifications of Eqs.(||) performing a 
small mass - m\ r and m 2r — expansion of the thresh- 
old functions. Using their asymptotic behavior one finds 
r\ = 0(u\ r , u 2r ) and : 

d t u lr = -euir + \(N + 8)uj r + 6u lr u 2r ] (6) 



d t u 2r = -eu 2r 



53 [* 



2r 



\2u\ r U 2r 



(7) 



These equations are those obtained directly in therjveak- 
coupling expansion of the TV- vector cubic model§0. For 
N =/= they display three nontrivial fixed points : 
the O(N) fixed point (u* lr > 0,u* 2r = 0), the TV- 
decoupled Ising one {u\ r = 0, u 2r > 0) and the cu- 
bic one (u\ r > 0,u 2r 7^ 0). The cubic fixed point 
is stable for N > N c and unstable for N < N c with 
N c (d = 4— e) = 4 — 2*|4- 0(e 2 ), as known from early per- 
turbative approachesElEl. For N — the only fixed point 
of the equations (||) and (0) is the Gaussian one but this 
is an artefact of the leading order of perturbation theory 
that disappears at secflnd order in e. This leads to the 
famous -y/e-expansionEllj. However, when the Eqs.(||), ob- 
tained with the simplest ansatz (|]) , are analyzed without 
e-expansion, the random fixed point already appears. All 
the fixed points found here persist when the dimension is 
lowered from d — 4 to the vicinity of d = 2, a case that 
deserves a special treatmentEj. 

Around d = 2, one also expects the Eqs.(@) to well 
describe, even quantitatively, the physics in the vicinity 
of the 0(N > 1) fixed-point. The reason is that, around 



this dimension, the phase transition associated to this 
fixed point takes place at a temperature of order e = d—2. 
Since the anomalous dimension is of the same order of 
magnitude, this situation must be well described even 
with our wild truncation that neglects high order deriva- 
tive terms. Around d — 2, l/(Nn r ) = T r provides the ef- 
fective perturbative temperature of the associated O(N) 
non-linear sigma (NLc) model. One also has to consider 
u 2 as a small perturbation around the O(N) fixed point. 
Thus, the procedure now consists in a double expansion 
of Eqs.(|J) in the large mass m\ r and in the small cou- 
pling constant u 2r . Using the asymptotic behaviors of the 
threshold functions, one obtains : 77 = T r /(2nN) and : 



d t T r 

dtUi r 

d t u 2r 



(N-2) 
C r 2n K 

-2u lr + 4u 2 lr ^ N ~ ^ In 2 



-2 + e H T r u 2r . 



(8) 
(9) 
(10) 



These equations admit a nontrivial fixed point with 
T* = 27re/(A-2), u* lr = n/(2(N- l)ln2) and u* 2r = 
which corresponds to the standard O(N) fixed point of 
the NLer model. It is then interesting to consider the flow 
of u 2r around this fixed point. One has : 



d t Su 2r 



(N + 6) 
(N-2) 



5u 2r 



(11) 



This shows that the cubic perturbation is relevant for 
N > N c and irrelevant for N < N c with N c (d = 2 + 
e) = 2 + 4e + 0(e 2 ), in agreement with the, perturbative 
approach of the NLcr model around d = 2c3. Concerning 
the fixed points with discrete symmetry (Ising O(l), iV- 
decoupled Ising, cubic and random Ising), they cannot, 
around d = 2, be properly described by our truncation. 
Actually, the large values of 77 at these fixed points around 
d = 2 likely require the inclusion in (^) of all powers of 
the field and high order derivative terms. 
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In d = 3, in order to obtain satisfying quantitative 
results, we have pushed the calculation using a more 
sophisticated ansatz than (^). We have considered all 
monomials up to order eight in the fields as well as all 
derivative terms including two derivatives and four fields 
and established the RG equations for the correspond- 
ing coupling constants. These equations are too long to 
be given hereEj and we just give our results. First, we 
have tested our calculation on the O(N) model. In the 
Ising and Heisenberg cases, we have found v — 0.627 and 
v = 0.719, respectively, whiieh are very reasonable com- 
pared to previous results (seea). The values of the anoma- 
lous dimensions r\ — 0.053 and r\ — 0.0452 are less convin- 
cing and related to the omission of higher order derivative 
terms in Tk- Considering the cubic anisotropy, we have 
determined an estimate of N c above which the cubic fixed 
point is stable and found N c (d = 3) = 3.1. This value 
differs significantly from those obtainedJxam 1 jf r l prl l t np r .- 
turbative approaches for which N c < ^Mt^B&B, 
For instance, in a receot. six- loop calculation, N c is found 
to be equal to 2.89(4)Ej, .a-isilue already obtained from 
lower order computationsESEZl. If this last result is con- 



firmed, our result should be imputed to the poor determi- 
nation of 77. The strong dependence of N c with respect to 
•q already appears in a hierarchical approximatiop_.where, 
neglecting any field renormalization, N c — 2.219E3. Note 
also that the result N c =p-|3.4 was found within the work 
of Newman and RiedelEj. For the RI fixed point we 
have found v = 0.67 and 7/ = 0.05 that give, through 
scaling laws, 7 = 1.306, (3 — 0.351 and a — —0.01. 
These rp L c jJ]| t j L ^7" J 1 !. £. nnrl agreement with perturbative 
methodsoBllj'Lj'tSLZrEj 7 and also with numerical esti- 
mates, while they show a certain discrepancy with the 
SF method that provides v = 0.697, 7 = 1.39 and 

7] = 0.015E3. 

Finally, we indicate that we have noted a rapid 
convergence of our results with respect to the order of 
the truncation in the fields. However, a challenging 
and still open question is that of the convergence when 
including higher order derivative terms which is essential 
for an improvement of the anomalous dimension. 

We thank M. L. Rosinberg, B. N. Shalaev and G. Tar- 
jus for discussions. 
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